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Abstract 

We study the depinning of a flux line by analytical and numerical methods 
applied to a phenomenological equation of motion. Transverse fluctuations 
do not influence the critical behavior of the longitudinal component, justify- 
ing "planar approximations". In an isotropic medium, longitudinal fluctua- 
tions have a roughness exponent = 1, and relax with a dynamic exponent 
z\\ K, 4/3; transverse fluctuations are suppressed {C,± = 1/2 < C||)i and relax 
more slowly, with z±_ = z\\ + 1. Anisotropy in the depinning threshold, or 
orientational dependence of force- force correlations, lead to new universality 
classes. 
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The pinning of flux lines (FLs) in Type-II superconductors is of fundamental importance 
to many technological applications that require large critical currents Upon applica- 
tion of an external current density J, the motion of FLs due to the Lorentz force causes 
undesirable dissipation of supercurrents. Major increases in the critical current density Jc 
of a sample are achieved when the FLs are pinned to impurities. There are many recent 
studies, both experimental and theoretical on collective pinning of FL's to point 
or columnar defects. Another consequence of impurities is the strongly nonlinear behavior 
of the current slightly above the depinning threshold, as the FLs start to move across the 
sample. Recent numerical simulations have concentrated on the low temperature behavior 
of a single FL near depinning mostly ignoring fluctuations transverse to the plane 

deflned by the magnetic fleld and the Lorentz force. Common signatures of the depinning 
transition from J < to J > Jc include a broad band (/^" type) voltage noise spectrum, 
and self-similar fluctuations of the FL proflle. 

The FL provides yet another example of a broad class of systems that exhibit a depin- 
ning transition; a dynamical critical phenomenon seen in Charge- Density Waves (CDWs) 
P^T^, interfaces in random media |]T3|,|T^, and contact lines There has recently 

been much theoretical progress in calculating the critical exponents of such transitions by 
renormalization-group (RG) methods. In this Letter, we extend these methods to study the 
full three-dimensional dynamics of a single FL at low temperatures. We show that fluctu- 
ations transverse to the average motion of the FL do not change the scaling properties of 
the longitudinal component. In turn, the transverse fluctuations are suppressed, and relax 
much more slowly (with a larger dynamic exponent). Thus, we not only formally justify 
ignoring transverse fluctuations (i.e. planar approximations) in numerical simulations near 
the depinning transition, but also demonstrate the rich anisotropic scaling properties of this 
system. 

In the simplest case of an isotropic pinning potential, the longitudinal fluctuations have 
roughness exponent C|| = 1, and dynamic exponent z\\ ~ 4/3, while C± = 1/2 < C|h 
z± = ^11 + 1 for the transverse components. When the pinning potential is anisotropic 
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and the FL is driven along a non-symmetric direction, transverse fluctuations are further 
suppressed to C± = 0, possibly with logarithmic deviations. If force-force correlations depend 
on the orientation of the FL as well as its position, the critical behavior falls into yet another 
universality class, possibly controlled by directed percolation clusters. Thus, the depinning 
behavior is quite rich, described by a number of distinct universality classes. 

We describe the shape of the FL at a given time t by r(x, t), where x is along the magnetic 
field B, and the unit vector ey is along the Lorentz force F. (See Fig. |l].) Point impurities 
are modeled by a random potential V{x,t), with zero mean and short-range correlations. 
In the presence of impurities and a bulk Lorentz force F, the energy of a FL with small 
fluctuations is, 

n = Jdx [lid^.rf + V (x, r(x, t)) - r(a;, t) ■ f} . (1) 

The simplest possible Langevin equation for the FL, consistent with local, dissipative dy- 
namics, is 

/^-^^ = -^ = 9lr + f {x, r{x, t)) + F, (2) 

where /i is the mobility of the FL, and f = —VrV. The potential V{x, r) need not be 
isotropic. For example, in a single crystal of ceramic superconductors with the field along 
the oxide planes, it will be easier to move the FL along the planes. This leads to a pinning 
threshold that depends on the orientation of the force. Anisotropy also modifies the line 
tension, and the elastic term in Eq. is in general multiplied by a non-diagonal matrix 
Kai3- The random force f (x, r), can be taken to have zero mean with correlations 

{fa{x, r)/^(x', r')) =6{x- x')Aayir - r ). (3) 

We shall focus mostly on the isotropic case, with AQ,^(r — r') = 5a^A(|r — r'|), where A is a 
function that decays rapidly for large values of its argument. 

While the flux line is pinned by impurities when F < Fc, for F slightly above threshold, 
we expect the average velocity v = |v| to scale as 
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vr^iF- F,f, (4) 

where P is the velocity exponent. Superposed on the steady advance of the FL are rapid 
"jumps" as portions of the hue depin from strong pinning centers. Such jumps are similar 
to avalanches in other slowly driven systems, and have a power-law distribution in size, cut 
off at a characteristic correlation length ^. On approaching the threshold, C, diverges as 

(5) 

defining a correlation length exponent z/. At length scales up to ^, the correlated fluctuations 
satisfy the dynamic scaling form, 

([r||(a;,t)-r||(0,0)]2) = |a;p^ii(?ll(t/|a;rii), (6a) 
t) - rx(0, 0)f) = \x\^'^^g^it/\x\'^), (6b) 

where (a and Za are the roughness and dynamic exponents, respectively. The scaling func- 
tions ga go to a constant as their arguments approach 0. Beyond the length scale ^, different 
regions of the FL depin more or less independently and the system crosses over to a moving 
state, described by different exponents C+,and z+. 

The major difference of this model from previously studied ones is that the position of 
the flux line, r{x,t), is now a 2-dimensional vector instead of a scalar; fluctuating along 



both ey and directions. One consequence is that a "no passing" rule [12|, applicable to 
CDWs and interfaces, does not apply to FLs. It is possible to have coexistence of moving 
and stationary FLs in particular realizations of the random potential. 

The effects of transverse fluctuations r± for large driving forces, when the impurities act 



as white noise, were studied earlier |T^, indicating a rich dynamical phase diagram. How 
do these transverse fluctuations scale near the depinning transition, and do they in turn 
influence the critical dynamics of longitudinal fluctuations near threshold? The answer to 
the second question can be obtained by the following qualitative argument: Consider Eq. (^ 
for a particular realization of randomness i{x, r). Assuming that portions of the FL always 
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move in the forward direction [0, there is a unique point rj_{x, ry) that is visited by the hne 
for given coordinates (x,r|j). We construct a new force field /' on a two dimensional space 
(x, r||) through f'{x,r\\) = /y (^x,r\\,r±{x,r\\)^. It is then clear that the dynamics of the 
longitudinal component r|j(a;,t) in a given force field f(a;, r) is identical to the dynamics of 
r||(x, t) in a force field f'{x, r||), with r± set to zero. It is quite plausible that, after averaging 
over all f, the correlations in /' will also be short-ranged, albeit different from those of 
f . Thus, the scaling of longitudinal fluctuations of the depinning FL will not change upon 
taking into account transverse fluctuations. However, the question of how these transverse 
fluctuations scale still remains. 

Certain statistical symmetries of the system restrict the form of response and correlation 
functions. For example, Eq. has statistical space- and time-translational invariance, 
which enables us to work in Fourier space, i.e. {x,t) {<1-,^)- For an isotropic medium, 
F and v are parallel to each other, i.e., v(F) = f (F)F, where F is the unit vector along 
F. Furthermore, all expectation values involving odd powers of a transverse component are 
identically zero due to the statistical invariance under the transformation r_|_ — t_i_. Thus, 
linear response and two-point correlation functions are diagonal. The introduced critical 
exponents are then related through scaling identities. These can be derived from the linear 
response to an infinitesimal external force field tu), 

Xcp{q,^) = (^~^~^ = ^c,f^Xa{q,uj), (7) 

in the {q,uj) — ^ (0,0) limit. Eq.(|D is statistically invariant under the transformation F — >■ 
F + e{q), r(g, uj) r(g, u) + q~'^e{q). Thus, the static linear response has the form X\\{qy ^ = 
0) = C(j = 0) = q^^. Since e\\ scales like the applied force, the form of the linear 

response at the correlation length ^ gives the exponent identity 

Cll + l/^ = 2. (8) 

Considering the transverse linear response seems to imply = However, the static 
part of the transverse linear response is irrelevant at the critical RG fixed point, since 



z\_ > z\\, as shown below. When a slowly varying uniform external force £{t) is applied, the 
FL responds as if the instantaneus external force F + e is a constant, acquring an average 
velocity, 

{dtr^) = Vo^iF + e)^ v^{F) + ^e,. (9) 
Substituting dv\\/dF^\ = dv/dF and dvj_/dFj_ = v/F, and Fourier transforming, gives 

Combining these with the static response, we see that the characteristic relaxation times of 
fluctuations with wavelength ^ are 
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(12) 



r^{q = r')-[<i'^) ~e^+^/^~r^ (13) 

which, using Eq.(^), yield the scaling relations 

/3=(^||-C||)^, (14) 
z± = + l/v. (15) 

We already see that the dynamic relaxation of transverse fluctuations is much slower than 
longitudinal ones. All critical exponents can be calculated from Cy, C^, and z\\, by using 
Eqs.d), (0), and (^). 

Equation (0) can be analyzed using the formalism of Martin, Siggia, and Rose (MSR) 
|[T9| . Ignoring transverse fluctuations, and generalizing to d dimensional internal coordinates 
X G 3?'^, leads to an interface depinning model which was studied by Nattermann, Stepanow, 
Tang, and Leschhorn (NSTL) and by Narayan and Fisher (NF) |]14[. The RG treatment 
indicates that impurity disorder becomes relevant for dimensions d < A, and the critical 
exponents in = 4 — e dimensions are given to one-loop order as C = e/3, z = 2 — 2e/9. 



NSTL obtained this result by directly averaging the MSR generating functional and 
calculating the renormalization of the force-force correlation function A(r), perturbatively 
around the freely moving interface [A(r) = 0]. NF, on the other hand, used a perturbative 
expansion of Z, around a saddle point corresponding to a mean-field approximation to Eq.(^ 
which involved tempore/ force-force correlations C{vt). They argue that a conventional 



low-frequency analysis is not sufficient to determine critical exponents. They also suggest 
that the roughness exponent is equal to e/3 to all orders in perturbation theory. 

Following the approach of NF, we employ a perturbative expansion of the disorder- 
averaged MSR partition function around a mean-field solution for cusped impurity potentials 



1^ . All terms in the expansion involving longitudinal fiuctuations are identical to the 
interface case, thus we obtain the same critical exponents for longitudinal fluctuations, i.e., 
C|| = e/3, z\\ = 2 — 2e/9 + 0(e^). Furthermore, for isotropic potentials, the renormalization of 
transverse temporal force- force correlations C±{vt) yields a transverse roughness exponent 
= 5^11 /2 — 2, to all orders in perturbation theory. The details of this calculation will be 
given elsewhere. For the FL (e = 3), the critical exponents are then given by 

C|| = l, Z||^4/3, u=l, 

(16) 

P ~ 1/3, a = 1/2, zi_ ^ 7/3. 

To test the scaling forms and exponents predicted by Eqs.(0) and (^), we numerically 
integrated Eq.(|D, discretized in coordinates x and t. Free boundary conditions were used 
for system sizes of up to 2048, with a grid spacing Ax = 1 and a time step At = 0.02. Time 
averages were evaluated after the system reached steady state. Periodic boundary conditions 
gave similar results, but with larger finite size effects. Smaller grid sizes did not change the 
results considerably. The behavior of v{F) seems to fit the scaling form of Eq.(P with an 
exponent /? ~ 0.3, but is also consistent with a logarithmic dependence on the reduced force, 
i.e., (3 = 0. The same behavior was observed by Dong et al. in a recent simulation of the 
1 + 1 dimensional geometry |^. Since z\\, and consequently (3, is known only to first order 
in e, higher order corrections are expected. By looking at equal time correlation functions 
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(see Fig. 0) , we find that transverse fluctuations are strongly suppressed, and that the 
roughness exponents are equal to our theoretical estimates within statistical accuracy. The 
excellent agreement for e = 3 suggests that the theoretical estimates are indeed exact. 

The anisotropy in critical exponents may be observed in a rectangular Hall geometry, by 
measuring the noise power spectra S\\ (u) and S_i{u!) of normal and Hall voltages, respectively. 
In a conventional type-II superconductor with point defects, at low temperatures and near 
depinning, Eqs.® suggest that P] Sa{uj) ~ cu""", where Oq, = {2(a + l)/2a — 1- Thus, 
Si\{uj) ~ uj~^^^^ whereas S±{uj) ~ 10^^"^ at small lu. 

The potential pinning the FL in a single superconducting crystal is likely to be highly 
anisotropic. For example, consider a magnetic field parallel to the copper oxide planes of a 
ceramic superconductor. The threshold force then depends on its orientation, with depinning 
easiest along the copper oxide planes. Eqs.(^, (1^) and ( pT]) have to be modified, since v 
and F are no longer parallel (except along the axes with r — r symmetry), and the linear 
response function is not diagonal. The RG analysis is more cumbersome: We find that, 
for depinning along a non-symmetric direction, the longitudinal exponents are not modified 
(in agreement with the argument presented earlier), while the transverse fiuctuations are 
further suppressed to (± = 2(\\ — 2 (equal to zero for (^y = 1) |2T|. Relaxation of transverse 



modes are still characterized by z± = z^^ + l/u, and the exponent identity (H) also holds. 
Surprisingly, the exponents for depinning along axes of refiection symmetry are the same as 
the isotropic case. 

Another feature of the problem is the possible dependence of the force-force correlations 
A on the orientation d^r of the FL in anisotropic superconductors. As in the case of interfaces 



22 1, such dependence leads to introduction of additional relevant terms in the MSR partition 



function, and invalidates the arguments leading to Eqs.(p|, p!^p!5D . The analogy to interfaces 
suggests that the longitudinal exponents for d = 1 are controlled by directed percolation 
clusters [^, with (\\ ^ 0.63. Since no perturbative fixed point is present in this case, it is 
not clear how to explore the behavior of transverse fiuctuations systematically. 

In conclusion, we studied the dynamical critical behavior of a single depinning FL in 
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a type-II superconductor at low temperatures. Using powerful symmetry arguments, we 
demonstrated the anisotropy in both the configurational and relaxational properties of the 
FL, which was confirmed by a formal RG treatment and numerical simulations. This justifies 
the "planar" approximation, widely used in numerical simulations of a depinning FL. Due to 
possible anisotropies in the pinning force, the depinning behavior is quite rich, encompassing 
a number of different universality classes. In the fully isotropic case, RG calculations suggest 
anisotropic scaling with z± = + 1, and (± = 1/2; the latter is confirmed by numerical 
simulations. Anisotropic potentials, and orientation dependent force correlations lead to a 
new exponents, some of which we have obtained. As the current density is further increased, 
other nonlinear terms become relevant and change the scaling of fluctuations. The effects of 
such terms were considered previously in Ref. fl^ for a single FL, and in Ref. p4[ for a dense 
collection of FLs. The critical behavior of an interacting ensemble of FLs near depinning is 
likely to be further influenced by complicated entanglement effects. 

We have benefited from discussions with O. Narayan and L.-H. Tang. This research 
was supported by grants from the NSF (DMR-93-03667 and PYI/DMR-89-58061), and the 
MIT/INTEVEP colloborative program. 
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FIGURES 



FIG. 1. Geometry of the fluctuating Flux Line. 

FIG. 2. Equal time correlations as functions of separation x, for a system of size 2048 at 
— Fc)/Fc ~ 0.01. The observed roughness exponents are close to the theoretical predictions of 
= 1, C± = 0-5, shown by solid lines for comparison. 
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